Quantum-interference-controlled three-terminal molecular transistors based on a 
single ring-shaped-molecule connected to graphene nanoribbon electrodes 
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We study all-carbon-hydrogen molecular transistors where zigzag graphene nanoribbons play the 
role of three metallic electrodes connected to a ring-shaped 18-annulene molecule. Using the nonequi- 
librium Green function formalism combined with density functional theory, recently extended to 
multiterminal devices, we show that the proposed nanostructures exhibit exponentially small trans- 
mission when the source and drain electrodes are attached in a configuration that ensures destructive 
interference of electron paths around the ring. The third electrode, functioning either as an attached 
infinite-impedance voltage probe or as an "air-bridge" top gate covering half of molecular ring, intro- 
duces dephasing that brings the transistor into the 'on' state with its transmission in the latter case 
approaching the maximum limit for a single conducting channel device. The current through the 
latter device can also be controlled in the far-from-equilibrium regime by applying a gate voltage. 

PACS numbers: 85.35.Ds, 85.65.+h, 73.63.Rt, 72.80.Vp 



Introduction. — The recent fabrication [TJ [2] of 
graphene nanoribbons (GNRs) with ultrasmooth edges 
has opened unforeseen avenues for nanoelectronics [3] by 
providing a novel type of semiconducting channel (GNRs 
with armchair edges of sub-10-nm width and band gaps 
> 0.4 eV) for field-effect transistors (FETs) pQ or low- 
dissipation interconnects [4] (metallic GNRs with zigzag 
edges). Furthermore, metallic GNRs have the capabil- 
ity to resolve one of the key challenges [5 for molecular 
electronics [6 — a well-defined molecule- electrode contact 
with high transparency, strong directionality and repro- 
ducibility. 

This is due to the fact that strong molecule-GNR 
7T-7T coupling makes possible formation of a continuous 
7r-bonded network across GNR and orbit als of conju- 
gated organic molecules [7]. In addition, unusual elec- 
tronic structure [4] of zigzag GNRs (ZGNRs) and the 
possibility to control it via edge doping [8 offer prospect 
for device designs that are outside of the scope of tra- 
ditional approaches in molecular electronics where single 
molecules are in contact with metals [5j |9] , such as gold 
or platinum which suffer from the lack of directional co- 
valent bonding thereby leading to poor reproducibility of 
most metal- molecule- metal junctions [10] , or carbon nan- 
otubes [11] where edges are not present and which do not 
have planar structure appropriate for aligning and pat- 
terning through lithographic methods suitable for high- 
volume production. 

Unlike carbon nanotubes (CNTs), which were em- 
ployed experimentally [11 as electrodes of molecular de- 
vices and whose contact to conjugated molecules has been 
investigated via first-principles quantum transport meth- 
ods [7] , surprisingly little is known about linear response 



or far-from-equilibrium transport properties of multiter- 
minal molecule-GNR hetero junctions. Although pristine 
ZGNRs have a band gap at low temperature due to edge 
magnetic ordering p~2j [13] , this gap is easily destroyed at 
room temperature [12 or by defects and impurities along 
the edge USE]. 

To obtain transistor effect, current through the molec- 
ular junction has to be controlled by the third elec- 
trode acting either as an electrostatically coupled gate [3 
in traditional FET designs [9] or as being attached to 
the molecule in quantum-interference-controlled transis- 
tor (QICT) concepts [15]. In the recently proposed 
QICT, source and drain electrodes are connected to a 
ring-shaped molecule in a configuration that ensures de- 
structive interference [16 of electron paths around the 
ring and, therefore, 'off' state of QICT with perfect zero 
leakage current. The third electrode, acting as the so- 
called Biittiker voltage probe which does not draw any 
current [16] , is then attached to the molecule to introduce 
dephasing and thereby 'on' state of QICT. 

The QICT also offers a playground for basic research 
on quantum interference effects [16] in molecular elec- 
tronics which, unlike in conventional two-dimensional 
electron gases or graphene, can manifest itself even at 
room temperature since molecular vibrations with de- 
phasing effect (i.e., those that change the length dif- 
ference between paths around the ring) can be sup- 
pressed [15] below 500 K. Such interference effects in 
transport through aromatic molecules have attracted 
considerable attention in the very recent experimental 
and theoretical studies [17] . Furthermore, potential ap- 
plications of QICT involve operation with greatly re- 
duced heat dissipation because the current flow is not 
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FIG. 1: (Color online) Schematic view of the pro- 
posed ZGNR|18-annulene|ZGNR three-terminal heterojuc- 
tions. The contact between the source and drain 8-ZGNR 
metallic electrodes and a ring-shaped 18-annulene molecule 
is made via 5-membered rings of carbon atoms (dark blue), 
while the electrodes are attached in a destructive configura- 
tion (1,9) for quantum interference. The source-drain current 
is controlled by a third electrode which is coupled to the de- 
vice either through bonding via 5-membered ring (a) or as 
an "air-bridge" top gate covering upper half of the molecu- 
lar ring (b). The hydrogen atoms (light yellow) are included 
to passivate the edge carbon atoms. The inset shows carbon 
atom numbering within the molecular ring, as well as three 
possible Feynman paths of a single electron entering the ring. 
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FIG. 2: (Color online) The effective source-drain transmis- 
sion T21 (E) for three-terminal QICTs shown in Fig. [l] and 
the corresponding LDOS at the Fermi energy. For compari- 
son, we also plot transmission T2i(E) for two-terminal devices 
in configurations (1,9) and (1,10) exhibiting destructive and 
constructive quantum interference, respectively. 



blocked by an energy barrier as in traditional FETs where 
it must be raised and lowered with each switching cycle 
in either silicon-based or carbon-based FETs [3]. 

Understanding of realistic junctions of ZGNR and 
molecules requires quantum transport methods combined 
with the first-principles input about atomistic and elec- 
tronic structure to capture charge transfer in equilibrium 
(which is indispensable to obtain correct linear response 
conductance of, e.g., carbon- hydrogen systems [18]) or 
charge redistribution (whose knowledge ensures the gauge 
invariance of the current-voltage characteristics) in the 
far-from-equilibrium regime driven by finite bias volt- 
age. The state-of-the-art approach that can capture all 
of these effects, as long as the coupling between the 
molecule and the electrodes is strong enough to dimin- 
ish Coulomb blockade effects [6], is the nonequilibrium 
Green function formalism combined with density func- 
tional theory (NEGF-DFT) [1 [18]. However, due to 
the lack of NEGF-DFT algorithms for multiterminal de- 
vices, previous efforts to model QICT based on benzene- 
gold [15] [19], annulene-gold [19] , benzene-CNT [19] and 
annulene-CNT [19] junctions have employed either semi- 
empirical models [15] (with truncated basis set and the 
connection between the leads and the molecule treated 
phenomenologically, thereby ignoring possible structural 
relaxation or hybridization between molecular tt states 
and contact a states) or NEGF-DFT [19] but in the ab- 
sence of third attached electrode. 

In this Letter, we analyze linear response 
and nonequilibrium transport properties of 
ZGNR| 18-annulene I ZGNR junctions, whose third 
electrode is made of ZGNR as well, by employing the 
very recently developed NEGF-DFT formalism for 



multiterminal nanostructures (MT-NEGF-DFT) [20]. 
The third electrode introduces dephasing by being 
attached to the molecule to act as an infinite-impedance 
voltage probe [16j, as shown in Fig. [TJa). Moreover, we 
find an even larger disruption of destructive quantum 
interference in a setup proposed in Fig. [IJb) where the 
third electrode plays the role of a top gate covering the 
upper-half of annulene ring while being separated from 
the two terminal device underneath by an air gap. Our 
principal result is shown in Fig. [2] where exponentially 
small transmission T2i(Ep) ~ 10 -6 between the source 
electrode 1 and the drain electrode 2 at the Fermi 
energy Ep in the junction with destructive geometry 
for connecting the two ZGNR electrodes increases by 
three orders of magnitude T2i(Ep) — 10 -3 with the 
attachment of an infinite-impedance voltage probe. 
Moreover, the junction becomes highly transparent 
T2i(Ef) — 0.4 due to induced states (Fig. |2| and charge 
transfer (Fig. |3| underneath the "air-bridge" top gate. 
We also analyze in Fig. [4] the gate voltage modulation 
of source-drain current in the far-from-equilibrium 
transport regime of device in Fig. []Jb). 

Device setup and MT-NEGF-DFT formalism.— 
Although naively one would expect that formation 
of continuous 7r-bonded network between carbon-based 
electrodes and conjugated molecules would ensure 
high contact transparency, early experiments [TT] on 
CNT I molecule I CNT hetero junctions have measured sur- 
prisingly small conductances for a variety of sandwiched 
molecules suggesting poor contact transparency. The 
first-principles analysis of different setups reveals that 
this is due to significant twisting forces when molecule is 
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connected to CNT via, e.g., 6-membered rings [7]. There- 
fore, to keep nearly parallel and in-plane configuration 
(hydrogen atoms of annulene slightly deviate from the 
molecular plane) of our ZGNR|18-annulene|ZGNR junc- 
tion, we use a 5-membered ring [7] in Fig. [I] to chemically 
bond ZGNR to annulene. The atomic structure of the 
hetero junctions in Fig.[l]is fully optimized by minimizing 
the atomic forces on individual atoms to be smaller than 
0.05 eV/A. This preserves the overall conjugation and 
leads to strong coupling and high contact transparency, as 
confirmed by Fig. [| where T 2 i(E F ) ~ 0.9 for ZGNR elec- 
trodes attached to annulene atoms (1,10) as an example 
of configuration with constructive quantum interference. 

The high contact transparency also makes it possible to 
apply the NEGF-DFT framework, which would otherwise 
be rendered insufficient due to electronic correlations [2T] 
emerging in the weak coupling regime that are beyond 
the mean-field DFT treatment [6 . In the NEGF-DFT 
formalism, the Hamiltonian is not known in advance and 
has to be computed by finding converged charge redistri- 
bution via the self-consistent DFT loop for the density 
matrix p = -jj^ J dE G < (E) whose trace gives the charge 
density \TE\ . The NEGF formalism for steady-state 
transport operates with two central quantities, retarded 
G(E) and lesser G < (E) Green functions, which describe 
the density of available quantum states and how electrons 
occupy those states, respectively. The technical details 
of the construction of the nonequilibrium density matrix 
p for multiterminal devices are discussed in Ref. [20] . 
In Fig. [3j we show the equilibrium charge density and 
the Hartree potential (obtained by solving the Poisson 
equation in the scattering region with the boundary con- 
ditions that match the electrostatic potentials of all three 
electrodes) computed through self-consistent loop for two 
QICTs from Fig. [j] 

In the coherent transport regime (i.e., in the absence 
of electron-phonon or electron-electron dephasing pro- 
cesses), the NEGF post-processing of the result of the 
DFT loop expresses the current flowing into terminal a 
of the device as: 
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Here the transmission coefficients Tp a (E,V a ,Vp) = 
Tr{l>(£, Vp)GT a (E,V a )G^} are integrated over the 
energy window defined by the difference of the Fermi 
functions f a (E) = {1 + exp[(E - E F - eV^/ksT]}- 1 
of macroscopic reservoirs into which semi- 
infinite ideal leads terminate. The matrices 
T a (E, V a ) = i[V a (E, V a ) - Et (E, V a )} account for 
the level broadening due to the coupling to the leads 
[E a (E,V a ) are the self-energies introduced by the leads 
whose electronic structure is assumed to be rigidly 
shifted by the applied voltage eV a ], thereby determining 
escape rates for electrons to leave the device and enter 
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FIG. 3: (Color online) The self-consistent equilibrium charge 
density [(a) and (b)] and the Hartree potential [(c) and (d)] 
computed within the central region of size of size 62 A x 
50 A of the three-terminal ZGNR|18-annulene|ZGNR devices 
shown in Fig.[l] The third electrode is either directly attached 
to the 18-annulene molecule [(a) and (c)], or it acts as "air- 
bridge" top gate at a distance 3.2 A away from the upper half 
of the molecular ring which it covers [(b) and (d)]. 



reservoirs where they are thermalized while the memory 
of their phase is lost. 

Our MT-NEGF-DFT code utilizes ultrasoft pseu- 
dopotentials and Perdew-Burke-Ernzerhof exchange- 
correlation functional. The localized basis set for DFT 
calculations is constructed from atom-centered orbitals 
(six per C atom and four per H atom) that are optimized 
variationally for the leads and the central molecule sep- 
arately while their electronic structure is obtained con- 
currently. 

Zero source- drain bias regime. — The interference ef- 
fects on quantum transport in traditional metallic and 
semiconducting mesoscopic rings are typically studied 
by applying magnetic flux through the ring to change 
the quantum phase accumulated by an electron circu- 
lating around the ring [16]. However, this method is 
not applicable to explore quantum interference effects in 
molecular-size devices, since magnetic field required to 
change the flux through the ring (e.g., from zero to h/e) 
far exceeds those available experimentally. Instead, one 
has to use particular geometry for attaching the elec- 
trodes to generate constructive or destructive quantum 
interference. For example, in the absence of the third 
electrode in device setups of Fig. [I] 7r-electron enter- 
ing the molecule at the Fermi level has the wavelength 
k F 1 2d (d is the spacing between carbon atoms within 
the molecule), so that for two simplest Feynman paths 
of length lOd (upper half of the ring) and Sd (lower half 
of the ring), the phase difference is 2k F d = tt (the two 
Feynman paths are illustrated in the inset of Fig. [I]). To 
get the full transmission T 2 \{E F ) between the source and 
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drain electrodes requires a sum over all possible paths, 
including more complicated ones [16 , but they all can- 
cel leading to a node T(Ep) = 0. We note that this 
single particle explanation is insufficient to understand 
interference effects in very small ring-shaped molecules, 
such as benzene weakly coupled to gold electrodes, which 
are often dominated by Coulomb blockade effect so that 
degeneracies between many-body states of the isolated 
molecule have to be taken into account [2T] . 

The reference curve in Fig. [2] is the zero-bias T2\(E) 
for the two-terminal device, which is exponentially small 
at Ep, but still a non-zero quantity. Although ex- 
act transmission node T2i(Ep) = is predicted for 
annulene-gold junctions using a phenomenological de- 
scription [15 , the NEGF-DFT framework finds that 
tunneling of a-electrons through hybridized a-orbitals 
in the gold electrodes and 7r-orbitals in the molecule 
can wash out such transmission nodes for not too large 
molecules [19]. In devices presented here, the high 
contact transparency facilitates injection of evanescent 
states [3 H3] from ZGNRs into the HOMO-LUMO en- 
ergy gap of the isolated molecule. The overlap of such 
states in the middle of the molecule generates a trans- 
mission resonance near Ep, while destructive interference 
for (1,9) connection of ZGNR electrodes superimposes an 
antiresonance dip T2\(Ep) — 10 -6 onto the two-terminal 
curve in Fig. [2] 

When the infinite-impedance voltage probe 3 is at- 
tached to the upper half of 18-annulene, new incoherent 
Feynman paths emerge (such as 1 — > 3 — » 2 illustrated 
in the inset of Fig. [I]) along which electrons propagate 
into the macroscopic reservoir through electrode 3 where 
they are dephased before entering the drain electrode 2. 
Such paths are not canceled, but instead generate [16] 
incoherent contribution T™(E) to the effective zero-bias 
transmission T2\{E)\ 



f 2 i(E) = T 21 (E)+T£ C (E) = T 21 (E)- 



T 2 3(E)T 31 (E) 



T 31 (E)+T 2 3(E)' 
(2) 

Figure pi shows that effective T2\{E) increases by three 
orders of magnitude due to the attached voltage probe. 
Furthermore, when the third electrode is not bonded, 
but instead covers half of the ring as shown in Fig. [TJb), 
T2\(E) c± 0.4 approaches the limit of a unit transmission 
through a single conducting channel (opened at Ep of 
ZGNR electrodes), which has been one of the key goals 
in molecular electronics pursuits of a perfect link [5]. 

The role of the "air-bridge" gate electrode in increas- 
ing the transmission by five orders of magnitude can be 
understood by examining the spatial profile of the local 
density of states (LDOS) in Fig. [2] and charge density in 
Fig. [3^b). They show how the top gate strongly modi- 
fies LDOS to enable charge transfer in the region under- 
neath. These effects are due to the fact that the top gate 
is positioned at the distance 3.2 A away from the planar 
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FIG. 4: (Color online) Transfer characteristics showing drain 
current I2 versus the gate-source voltage V gs for the device 
in Fig. [TJb) biased by the drain-source voltage Vds — 0.2 V. 
Right axis gives values of current in the third (top gate) elec- 
trode. The inset plots I2 vs. Vds for the two-terminal device. 



device, so that hybridization due to overlap of gate or- 
bitals with atomic orbitals underneath generates LDOS 
directly connecting leads 1 and 2. 

Finite source- drain bias regime. — Driving mesoscopic 
and nanoscopic devices out of equilibrium by applying a 
finite bias voltage typically corrupts their quantum co- 
herence (e.g., electrons at different energies within a bias 
voltage window have different interference patterns), so 
that standard interference effects are observed in the lin- 
ear response regime [16] . Concurrently, at finite bias volt- 
age Vds between the source and drain electrodes, the at- 
tached third electrode as in the device in Fig. [TJa) would 
draw finite current [15], which is an example of leakage 
currents whose minimization is one of the major tasks in 
downscaling FETs in digital electronics [3]. Here we use 
the setup in Fig. [TJb) to investigate the range of mod- 
ulation of source-drain current at Vds — V2 — V\ — 0.2 V 
via the top gate positioned at the distance 5.3 A away 
from the device. Figure [4] shows that current can be 
modulated from very small to finite value by changing 
the gate voltage V gs , while the selected distance ensures 
that tunneling leakage current into the top electrode is 
very small. Thus, this type of a device could act as ei- 
ther a molecular-scale FET [9] or as a bipolar junction 
transistor capable of amplifying current in the third lead. 

Conclusions. — In conclusion, using NEGF-DFT for- 
malism, recently extended to multiterminal devices [20] , 
we analyzed quantum transport through ZGNR| 18- 
annulene I ZGNR junction in the presence of the third 
ZGNR electrode attached either as the infinite- 
impendence voltage probe or serving as the "air-bridge" 
top gate covering half of the molecular ring. While both 
choices for the third electrode disrupt destructive quan- 
tum interferences (responsible for exponentially small 
transmission in the corresponding two-terminal device), 
transmission increases substantially further when using 
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the "air-bridge" top gate which enhances the LDOS un- 
derneath to directly connect the source and the drain. 
Moreover, the latter device in its 'on' state even ap- 
proaches the long-sought limit [5] in molecular electronics 
of a unit transmission through an organic-based device 
with a single conducting channel, because of high con- 
tact transparency brought by the proposed application 
of graphene-based electrodes. 
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